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ABSTRACT 

Solar sails provide a propellant free form for space propulsion. These are large flat 
surfaces that generate thrust when they are impacted by light. When attached to a space vehicle, 
the thrust generated can propel the space vehicle to great distances at significant speeds. For 
optimal performance the sail must be kept from excessive vibration. Active control techniques can 
provide the best performance. However, they require an external power-source that may create 
significant parasitic mass to the solar sail. However, solar sails require low mass for optimal 
performance. Secondly, active control techniques typically require a good system model to 
ensure stability and performance. However, the accuracy of solar sail models validated on earth 
for a space environment is questionable. An alternative approach is passive vibration techniques. 
These do not require an external power supply, and do not destabilize the system. A third 
alternative is referred to as semi-active control. This approach tries to get the best of both active 
and passive control, while avoiding their pitfalls. In semi-active control, an active control law is 
designed for the system, and passive control techniques are used to implement it. As a result, no 
external power supply is needed so the system is not destabilize-able. Though it typically under- 
performs active control techniques, it has been shown to out-perform passive control approaches 
and can be unobtrusively installed on a solar sail boom. Motivated by this, the objective of this 
research is to study the suitability of a Piezoelectric (PZT) patch actuator/sensor based semi- 
active control system for the vibration suppression problem of solar sail booms. Accordingly, we 
develop a suitable mathematical and computer model for such studies and demonstrate the 
capabilities of the proposed approach with computer simulations. 

INTRODUCTION 

Solar sails can propel a space vehicle using thrust generated when energy of photons 
from the sun are reflected from large flat surfaces. Though the force generated is small, 
propellants are not required and the space-craft can achieve large velocities because the external 
force acts for the life of the mission ([1],[2]). A Solar Sail configuration to address NASA’s future 
propulsion needs is under development by a team from L’Guarde, Jet Propulsion Laboratories, 
Ball Aerospace, and NASA Langley Research Center ([1], [2]). This 100m square sail 
configuration consists of four identical quadrants each containing a triangular sail. The four sails 
(made from Mylar) are supported by four booms (at right-angles to each other) attached to the 
solar sail craft. Control vanes (small reflective areas) that resemble a scaled quadrant are 
attached to the free tip of each boom via actuators. Thrust is vectored by rotating the vanes 
around the boom axis. Full three axis control is achieved by actuating three (or four) vanes 
simultaneously [2], 

The efficiency with which energy from the reflected photons is converted into solar sail thrust will 
be affected by sail flatness. Since solar sails have no inherent bending stiffness, and little mass, 
sail dynamics will be dominated by boom dynamics. In addition, the booms are long, flexible, and 



lightly damped structures. Therefore, disturbances to the boom (typically a truss) can result in 
large amplitude, long duration vibration motion that will have a debilitating impact on the 
performance [3], Minimizing the adverse effects of this vibration will require vibration mitigation 
technologies. Piezo-electric (PZT) actuation devices can convert mechanical vibration energy to 
electrical energy and vice versa. They also have wide frequency response spectrum and strong 
electro-mechanical coupling. If bonded to the surface of a flexible structure, they are conceptually 
capable of operation in the space environment and can therefore be classified as “space- 
realizable” [3], Furthermore, thin, light-weight and flexible PZT actuators such as polyvinyledine 
(PVDF) film and a macro-fiber composite patch can be packaged in an inobtrusive manner on (or 
within) the truss such as the L’guard Truss without introducing significant mass loading. For these 
reasons, PZT devices have been extensively used as sensors and actuators for active and 
passive vibration suppression of flexible structures [4], 

Active control of vibration involves utilization of powered actuation devices and easily benefits 
from modern control design techniques. Consequently it can potentially provide the best 
performance for a given hardware configuration. However, it requires a significant external energy 
source to drive the actuators and can destabilize the structure if improperly designed [5], 
Furthermore the required energy source will introduce significant parasitic masses, and thus this 
approach is not well suited to solar sail craft where low masses are critical. Passive approaches 
have the advantage that they can not destabilize the system under any circumstances because 
they only take vibration energy out of the structure. In this approach a Piezo-electric actuator is 
shunted (e.g. with a resistor or resistor and actuator). As a result, mechanical vibration energy is 
converted to electrical energy, which is then dissipated via a passive circuit. This approach 
typically under-performs active approaches and may require large inductance (i.e. a heavy coil) 
for low frequency vibration suppression. Since the boom is long and lightweight, the dominant 
modes will have a low frequency. Therefore a heavy coil may be required for each actuator. Since 
several actuators may be needed, the large total mass of all the coils makes passive control with 
PZT devices un-ideal for ultra-light weight structures like solar sail craft. 

Recently a new vibration suppression approach called “semi-active” or “hybrid” vibration 
suppression has been introduced. The idea is to convert vibration energy to electrical energy via 
PZT actuators attached to the vibrating structure, accumulate the electrical energy, and then 
dissipate the stored energy in a controlled manner. For example, a resistive, capacitive or 
inductive shunt circuit (for the PZT actuator) is switched on and off according to the phase of a 
vibrating mode of the structure [5], Since vibration energy is always extracted from the system, 
this approach cannot destabilize the system. A benefit of this approach is the fact that external 
power to drive the actuators is not needed. However, a small external energy source may be 
needed to power the switching circuit. In view of the above, semi-active control techniques are 
well suited to vibration suppression of a solar sail boom. However, research on this approach is 
limited and further studies are needed. 

In Reference 6, we showed that a solar sail boom could be modeled as an Euler-Bernoulli beam 
with clamped-free boundary conditions. As an extension of this effort, we herein use an Euler- 
Bernoulli beam to illustrate the capabilities of semi-active vibration suppression for solar sail 
booms. We begin by developing a mathematical model of a classical Euler-Bernoulli beam with 
clamped-free boundary conditions and integrated PZT actuators/sensors. The model can 
accommodate multiple PZT actuator/sensor patches at arbitrary locations on the beam. Next we 
validate the model with experimental input-output data. Since models of solar sail booms are 
expected to contain significant uncertainties, we consider an active control law (in this case 
velocity feedback) that is inherently robust to model uncertainties. Semi-active control is then 
simulated as a switching logic that causes the charge generated in the PZT actuators (due to 
beam vibration) to have the same sign as the charge that would be induced in the actuator by the 
active control law in a single degree of freedom system. The results of the simulation are 
provided and compared with equivalent results if a passive approach had been used. 


EXPERIMENTAL SETUP 



The test setup is shown in Figure 1. A galvanized steel beam was clamped at one end and free 
to move at the other end. A pair of Macro Fiber Composite (MFC) PZT patches (manufactured 
by smart material corporation) were glued to opposite sides of the beam close to the clamped end 
and wired so that their individual response to an identical voltage will be 180° out of phase. A 
laser displacement sensor was used to measure the beam displacement at the farthest point on 
the patch from the clamp (i.e. 56mm from the clamp). Additional details are provided below. 



Figure 1: Test Setup 


PROPERTY 

SYMBOL 

VALUE 

Modulus of Elasticity of the (Steel) Beam 

E b 

207 x 10 s N/M' 

Modulus of Elasticity of the PZT 

Ep 

30.336x 10* N/M^ 

Density of the (Steel) Beam 

Pb 

7850 Kg/M J 

Free Strain per volt (average of: Low field ~ High field) 

P 33 

0.75 ~ 0.9 ppmA/olt 

Apparent Internal Capacitance of the PZT Actuator 

Ca 

1.7052 x10' s F 

PZT Actuator Input Voltage Range 


-500V ~+1 500V 


Table 1: Material Properties 


PARAMETER 

SYMBOL 

VALUE 

Length of the beam (from clamped point to free end) 

Lb 

288mm 

Beam Width 

W b 

25.4mm 

Beam Thickness 


0.508mm 

Location of each MFC patch from clamped end 


16mm - 56mm 

Length of each MFC patch 


40mm 

Width of each MFC patch 


25.4mm 

Length of active area (Epoxy + PZT at center) of MFC patch 

Lp 


Width of active area (Epoxy + PZT at center) of MFC patch 

W P 

14mm 

Thickness of active area (Epoxy + PZT at center) of MFC patch 

tp 

29mm 

Thickness of epoxy alone (around the edges) of MSFC patch 


14mm 


Table 2: Dimensions 


TEST PROCEDURE 

To obtain validation data for the physics model, a 1 Hz square wave with amplitude 1 Volt was 
generated with a National Instruments data acquisition system and then amplified (with a Trek 
Model PZD 700-2-L-CE Piezo Driver/Amplifier) by a factor of 200. A laser displacement sensor 
(ILD 1800-200 made by Micro Epsilon) was used to measure the displacement of the Beam/PZT 
at 56mm from the clamp point at a sample rate of 1000Hz. The Frequency response functions of 
the input-output data were computed with the MATLAB spafdr function, while sixth order models 
were identified with the n4sid in the MATLAB System Identification toolbox. The validated model 
was subsequently used for the simulations presented herein. In this study the vibration 
suppression schemes will be designed to target the first vibration mode. Accordingly, an 

“optimal” inductance L for passive vibration suppression is I = 1/ co 2 C a where 0 ) is the natural 
frequency of the first mode in radians per second, and C a is the inherent capacitance of the PZT 
actuator. Similarly, an optimal value of inductance for semi-active control isi = 0.1 / co 2 C a . 


MATHEMATICAL MODEL 


Assumptions. 

1. Beam has a uniform mass distribution and a uniform stiffness, that are m(s) = m and 
EI(s ) = El respectively. 

2. Perfectly straight beam. 

3. External forces and moments do not induce any torsional effects. 

4. Cross-sectional dimensions are small compared to beam length. Note that shear and rotary 
inertia effects can be neglected under this assumption. However, these effects become significant 
as mode number increases. 
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5. Small bending deformations, specifically: 
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6. The beam is stiff in the axial direction, so length change due to axial forces is negligible. 


EQUATION OF BENDING MOTION 

We considered a beam clamped at one end and free at the other end and derived equations of 
transverse motion y(s,?)in response to applied forces and moments. At rest, the beam is 
vertical to minimize gravity effects and thereby simulate the zero gravity effects of space. Its 
system equation for bending motion is ([6] and [7]): 
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where s is the longitudinal spatial coordinate, y is the transverse displacement, El is the 
bending stiffness, T{L ) is axial tension at s = L , m is the mass per unit length, L is the length 

of the beam, t is time, w is the spatially distributed transverse force and M is the applied 
moment per unit length. 

We used admissible functions in our analysis because of difficulties computing the mode 
shape functions 0 . (s) . The admissible functions are the mode shape functions 

of the classical Euler-Bernoulli Beam, i.e: 

Uj (s) = c y sin (pjs) + c 2j cos(PjS) + c 3J sinh(/^) + c 4J cosh (fijs) (2) 


where £> 7 , and c Yj , c 2J , c 3j , c 4/ , j- 1,2,3,.- are computed from boundary conditions. 


BENDING MOMENTS DUE TO PZT PATCHES 


1. 



y Axis 


if]?ZT 


tp 

4 







T 


X Axis 


Beam 






T 


Figure 2: Cross-sectional view of PZT Actuator Pair on a Beam 


The active area of the PZT patches on our test-bed is not as wide as the beam even though the 
PZT patch and the beam have the same width. Though many texts (e.g. [8]) provide physics 
equations for a beam with integrated PZT actuators, we are not aware of any results that explicitly 
accounts for the case where PZT patches have a different width from the beam. Motivated by 







this, we developed physics models for the case where several such PZT patches are attached to 
arbitrary locations along an Euler-Bernoulli beam. Our models consider the case where the PZT 
patches are used as actuators or sensors. Such models capture the dual actuator/sensor nature 
of PZTs needed for our study. 

Figure 2 shows a cross-sectional view of a beam with thickness t b that has two PZT patches 
(each with thickness ) at opposite sides of the beams longitudinal axis. The strain s p in the 
PZT is: 

£ — ^33 

p T ' 

- ^ p 

where V A is the voltage applied to the MFC PZT patch across its 3 direction. We assume a 
linear strain distribution s(z) across the Beam/PZT cross-section induced by the PZT patch pair 
is given by 

s(z ) - a z . 

Note that stress is not continuous though strain is. From Hookes’ law: 

a b (z) = E b s(z) = E b a z 

where E b is the Elastic modulus of the beam, For the case where the voltage applied to the top 
PZT is 180° out of phase with that applied to the bottom PZT, we get: 


Top PZT: a p =E p a z-E p s p =E p (a z-s p ) 

Bottom PZT: a B p =E p a z + E p s p =E p (a z + s p ) 


where E p is the Elastic modulus of the PZT. For equilibrium about the center of the beam, we 
have 
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We can solve for a by integrating the expression above to get 
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Therefore the induced moment per unit length M in the beam due to the PZT patches is: 


M At) = E b I b a(t) = E b I b 
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is the moment of inertia of the beam. Next we calculate 


j&O) 




d 2 M p (s,t) 
ds 2 

which is needed to compute transverse forces induced in the beam by the PZT patches. Let 
M ■ (s, t ) denote the moment due to the j th PZT patch (from a set of m 2 actuator pairs) located 

ata ; < x<bj , and V A ( t ) is the voltage applied to the actuator, then 

77? 2 

M„(s, 0 = 

7=1 

with 

M b 0,0 = (r)(t/(s - a,) - t/(i -i y )) 

where is a constant implicitly defined in eq. (3) and £/(’) is the unit step function. From the 
definition of Dirac’s delta function, we have 

o 

U(s-z)= J 8(x-z)dx . 
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Accordingly for i=1,2,...,n : 
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VOLTAGE GENERATED BY PZT PATCHES 
From eq. (3), we have 


M D (t) = K VA V A (t). 


(4) 


It is well known that 



M = 


El 


where r is the radius of curvature of a beam in flexure. Also we can set 

1 _ d 2 y{x) 
r dx 2 

assuming small deformations. Applying these two expressions to our system, we have 

Equating the RHS of the expression above with the RHS of eq. (4) gives 

d 2 y{x ) 


E h L 


dx 2 


= K VA V A (t). 


Assuming uniform flexural deformation due to the i th PZT at a t <x< b ( , we have 

d 2 y(x) 


EJ, 


dx 1 


K VA V Ai (t) 


For non-uniform flexure, we integrate the expression above from a, to 6, and get 


E b I t )^^-dx= \K VA V M (x,t)dx. 


dx 

Uj Llj 

Next we assume the voltage-to-curvature equation above can be inverted as follows: 
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Where V G ( t ) is the voltage generated in the I th PZT (located at ct t < x < b t ) by flexure 

d 2 y(x,t ) 


dx 2 


acting within a, < x <b r Assuming separation of variables: 


y{x,t) = Y d $ j {x)q ] (f) 


we have 

v G ,{t) = 


E b I b 
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Note that the preceding expression can be used for an arbitrary number of PZT patch pairs each 
of which can have arbitrary locations, width and thickness. 
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SIMULINK MODELS 

We consider the case where the beam is initially excited by a 100 Volt external signal that acts on 
the transducer from time t=0.01 to t=0.01 1 seconds. Once excited, the resulting vibration is 
suppressed using either a passive or semi-active approach. This scenario was simulated in 
Simulink using Powergui (a graphical user interface for the analysis of circuits and systems). In 
this environment, the beam was represented by a transfer function computed from the physics 
model. For ease of exposition, we considered the collocated actuator/sensor case for which 



velocity feedback is an appropriate active control strategy. Onoda et al. [5] demonstrated that 
semi-active vibration suppression is achieved by periodically shunting the terminals of a PZT 
transducer actuator across a series circuit consisting of an inductor and resistor according to the 
following simple rule. First design an active control signal (e.g. velocity feedback) and use it to 
compute the desired control signal U Q . Next periodically discharge the charge Q built up in the 
PZT (due to vibration) such that Q U Q > 0 whenever U Q > 0. Accordingly, the following switching 
rule was devised to achieve this: 

T urn on the switch when U Q V < 0 and turn it off when U Q V > 0, 

where V is the voltage across the Transducer terminals. This switching logic was implemented in 
the Simulink diagram in Figure 2. Passive control was subsequently simulated by eliminating the 
switch in the Simulink diagram of Figure 2. 


| Continuous 
pwergui 



To Vfoitepace 



Figure 2: Simulink Model of the Beam under Semi-Active Control 

RESULTS AND DISCUSSION 

Onoda et. al [5] successfully implemented semi-active vibration control methods using PZT stack 
actuators. In their systems, the electrical energy generated when a PZT stack actuator is 
attached to a vibrating structure was dissipated passively in a manner that yields more vibration 
suppression than conventional passive methods. In our study, we plan to evaluate the 
effectiveness of their approach for structures with multiple integrated MFC PZT patches each of 
which has significantly lower inherent capacitance than conventional PZT stack actuators. As a 
first step, we developed a multi-actuator model and used it to study a beam with a co-located 
actuator pair. 



MODEL VALIDATION RESULTS 


First we consider the case where the PZT acts like a sensor. To study this, we used an impulse 
hammer to strike the beam at its free end and measured its response. The applied input force, 
transverse displacement, and generated voltage from both PZT actuator/sensors are provided in 
Figure 3. The fact that an impulse struck the beam around 2.2 seconds is evident from the PZT 
voltages in the bottom two plots of Figure 3, but not apparent from the Hammer and Laser 
Displacement plots at the top of Figure 3. Unfortunately, the magnitude of the applied force was 
apparently below the resolution of the impulse hammer and could therefore not be measured 
because of the noise. Similarly, a significant amount of ambient noise was present in the Laser 
displacement sensor signal. 
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Figure 3: Input-Output Data for Impulse Force 

The quality of the correlation between this set of input-output data was poor due to the presence 
of significant noise. This is evident from the incoherent appearance of the Frequency Response 
Function (FRF) computed from the raw data and the low magnitude and flat profile of the 
identified models. Consequently, the FRFs of Force-Input-to-Displacement-output of Figure 4, 
and the Force-input-to-Voltage-output of Figure 5 suffer from similar problems. However, the 
location of resonant peaks of the experimental data and the mathematical model appear to 
roughly match. 

Next we obtained data for the case where the PZT acts as an actuator. Unfortunately, we had 
the same ambient noise problems again. The resulting FRF depicting the relationship between 
the applied voltage and beam displacement is depicted in Figure 6. In all cases the physics 
model had first second and third vibration mode natural frequencies of 5.0408, 32.3952, 89.8178 
Hz respectively and a damping ratio of 0.01 was arbitrarily assumed for all three modes. 









Figure 4: Frequency Response Plots for Force Input and Voltage Output 



Figure 5: Frequency Response Plots for Force Input and Displacement Output 
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Figure 6: Frequency Response Plots for Voltage Input and Displacement Output 


VIBRATION SUPPRESSION RESULTS 

The results of the simulation of passive and semi-active vibration suppression are shown in 
Figure 7. The top part of Figure 7 shows the open loop (solid line) and closed loop (dashed line) 
response under passive control. Similarly, the bottom part of Figure 7 shows the open loop (solid 
line) and closed loop (dashed line) response under semi-active control. The superior 
performance of the semi-active control scheme is readily apparent. 
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SUMMARY AND CONCLUSIONS 


We demonstrated semi-active control with PZT actuators may be well suited to the vibration 
suppression problem of a long slender beam designed for solar sail craft. If a decentralized 
approach is used, there will be no need to run long cables and the system will be robust against 
failure of some of its sensors, actuators or controllers. Since only small external power is needed 
to power the control logic, the entire vibration control system should impose a much smaller 
parasitic mass on the truss than an equivalent active and passive approach. 

FUTURE WORK 

In the future, data with significantly less noise will be used to validate the physics model. 

Obtaining such data will require a more suitable test-bed. Also closed loop experimental data is 
needed to fully evaluate and confirm the suitability of MFC PZT patch actuators for the vibration 
suppression problem of solar sail booms. Finally, the issue of how to implement the shunting 
circuits in a space environment should be investigated. 
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